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ABSTRACT 



An iterative technique for the retrieval of temperatures 
at each of 100 levels ranging from 1000 mb to 0.01 mb is 
evaluated. Clear-column radiance data in the carbon dioxide 
channels of the VTPR of NOAA 2 are used in inverting the 
radiative transfer equations to deduce the T(P) profile. The 
retrieval technique includes the computation of atmospheric 
transmittances due to three atmospheric absorber masses 
(carbon dioxide, water vapor, and ozone) and non-homogeneous 
temperature-pressure effects along the vertical. The program 
also corrects these transmittances for zenith path differences 
between the satellite and the retrieval site when the site is 
not directly below the sensor. A normalized wave number, 
close to 700 cm \ but decreasing with decreasing values of 
pressure, was used in solving the Planck function for tempera- 
ture, level-by-level. The retrieved temperatures were then 
compared with matching radiosondes at 15 key atmospheric 
levels and the resulting errors analyzed on a globally averaged 
and latitude band basis at each of the 15 levels. 



4 



TABLE OF CONTENTS 

I. INTRODUCTION 12 

II. SATELLITE DATA 17 

A. NOAA 2 VTPR DESCRIPTION 17 

B. CLEAR-COLUMN PROCEDURE 19 

III. TRANSMITTANCE MODELS IN THE VTPR CHANNELS 26 

A. ATMOSPHERIC TRANSMITTANCE COMPUTATIONS: 

GENERAL CONSIDERATIONS 26 

B. WATER VAPOR - MIXING RATIO 

ESTIMATION PROCEDURE 28 

C. ADJUSTMENT OF WATER VAPOR AMOUNTS FOR 

NON-HOMOGENEOUS ATMOSPHERES 31 

D. WATER VAPOR TRANSMITTANCE MODEL 35 

E. CARBON DIOXIDE TRANSMITTANCE MODEL 38 

F. TRANSMISSIVITY ADJUSTMENT FOR ZENITH ANGLE - 40 

IV. RETRIEVAL OF TEMPERATURE PROFILES 43 

A. INPUT DATA USED FOR TEMPERATURE RETRIEVAL — 43 

1. First Guess Temperatures 43 

2. Other Input Data 4 5 

B. MATHEMATICAL DEVELOPMENT OF THE 

TEMPERATURE RETRIEVAL 45 

C. STEPWISE TEMPERATURE RETRIEVAL 52 

D. PROGRAM OUTPUT 55 

V. RESULTS OF THE VTPR TEMPERATURE RETRIEVAL 58 

A. GLOBALLY AVERAGED RESULTS 58 

B. LATITUDE BAND RESULTS 60 

VI. CONCLUSIONS 64 



5 



APPENDIX A-l, A-2 , . . . ,A-15. Histograms of the 

Globally Averaged Errors in Temperature, 

T(P) , as Retrieved at the Levels 
P = 1000, 850, 700, 500, 400, 300, 250, 

200, 150, 100, 70, 50, 30, 20, 10 mb 66 

APPENDIX B Ozone Transmittances for the Six 

VTPR Clear-Column Channels 81 

APPENDIX C U. S. Standard Atmosphere Supplements 84 

COMPUTER PROGRAM 86 

A. MAIN PROGRAM 86 

B. SUBROUTINES 100 

SAMPLE COMPUTER PROGRAM OUTPUTS 109 

A. CARBON DIOXIDE TRANSMITTANCES 109 

B. WATER VAPOR TRANSMITTANCES 112 

C. TOTAL ATMOSPHERIC TRANSMITTANCES 115 

D. SAMPLE OUTPUT FROM STATION 1 118 

LIST OF REFERENCES 119 

INITIAL DISTRIBUTION LIST 122 



6 



LIST OF TABLES 



1. VTPR channel designators and centroid wave numbers 
for filter set number four of NOAA 2 VTPR 

instrument number 1 (used in March 1973) 19 

2. Eleven levels used for predictor levels in water 

vapor retrieval 30 

3. 56 level atmosphere with merged climatology 

and 15 level first guess field for sample sounding — 44 

4. Normalized wave numbers and corresponding pressure 
levels used in inverting the Planck function 

for temperature retrieval 51 

5. Latitude bands used for error analysis of the 

VTPR retrieved temperatures 57 

6. Retrieved temperature error on a global basis 59 

7. Retrieved temperature error on a latitude 

band basis 61 



7 



LIST OF FIGURES 



1. Derivative of transmittance with respect to 

(pressure) 2/7 for the SIRS-B radiometer channels 13 

2. Satellite tracks for NOAA 2 VTPR coverage 18 

3. VTPR scan pattern and data analysis array 21 

4. NOAA 2 VTPR lateral scan geometry at 

raw radiance spots 22 

5. Scanning patterns used by the National Environmental 

Satellite Service (NESS) for pair comparison in 
determining VTPR clear-column radiance 24 

6 . Multi-level atmosphere showing P(J), T(J), U(J) 32 

7. Transmittance curves for atmospheric layers 1 and 2 

illustrating modification to effective U(J) in a 
non-homogeneous atmosphere 34 

8 . Thirty three K-layer atmosphere model used for the 

integration of the radiative transfer equation 47 

9. Carbon dioxide transmittances in channels 

1 through 6 50 

10. Weighting functions used in determining 

a weighted mean Planck function 51 



8 



TABLE OF SYMBOLS AND ABBREVIATIONS 



B i [T (P) ] 


Planck radiance function for temperature T 
at pressure level p and channel i 


B ± (K) 


Layer mean Planck function for layer K, in 
channel i 



/ n \ 

B' ' [v,T(J)] Weighted mean Planck function at iteration n 

and for level J, using the normalized wave number v 



CLRAD 


National Environmental Satellite Service (NESS) 
program for obtaining clear-column radiances 


ITPR 


Infrared Temperature Profile Radiometer 


J 


The subscript used to denote the levels in a 
100 level atmosphere where 

p^/^ _ (o 01) 

J = 1 + 99 [- 277] 

(1000) 7 ' - (0.01) 7 


K 


The subscript used to denote the layers in a 
33 layer atmosphere 


mb 


millibar 


NASA 


National Aeronautics and Space Administration 


NESS 


National Environmental Satellite Service 


NMC 


National Meteorological Center 


NOAA 


National Oceanic and Atmospheric Administration 


NTP 


Normal Temperature and Pressure conditions 


V i 


Wave number at the center of channel i 


V 


• —1 , 

Normalized wave number, v = 700 cm and is a 

decreasing function of pressure 


P(J) 


Pressure at level J 


P s 


Pressure at the surface of the earth (1000 mb) 


RTE 


Radiative transfer equation 



9 



SIRS 


Satellite Infrared Spectrometer 


SR 


Scanning Radiometer 


T [P ( J) ] 


Temperature of the atmosphere at the pressure 
level P(J) 


T ± (J) 


Atmospheric transmittance at level J, in channel i 


T [X ( J ) ] 


Atmospheric transmittance at level J from 
constituent X 


ym 


_6 

micrometers = 10 meters 


VTPR 


Vertical Temperature Profile Radiometer 


* 


Fortran symbol for multiplication used, if 
necessary, for clarification 



10 



ACKNOWLE DGEMENTS 



The author wishes to express his appreciation to Professor 
Frank L. Martin for his generous assistance, infinite patience 
and guidance in the research and development of this paper. 

Grateful acknowledgements are also extended to Dr. D. Q. 
Wark and Dr. H. E. Fleming of NESS for making available the 
archival tapes of the select data of March 1973, which made 
this paper possible. 

It is appropriate here to recognize the generous support 
of Dr. L. M. McMillin and Dr. M. P. Weinreb who provided not 
only operational programs which were consistent with the 
instrument in use, but also, when requested, explained the 
physical rationale of such programs in connection with 
retrieval procedures at NESS. 



11 



I. INTRODUCTION 



The quantitative measurement of temperatures at known 
pressure levels in the atmosphere is fundamental to present 
day weather analysis and forecasting. Until the launch of 
the Nimbus 3 satellite with its Satellite Infrared Spectrometer 
(SIRS) in April 1969, free air temperature soundings were 
obtained entirely from balloons, aircraft, and rockets. Since 
these sources were readily available only from populous land 
areas and a few ships, large areas of the world were without 
any reports. World-wide distribution of temperature sounding 
data was made possible by the use of satellites. The retrieval 
of a temperature profile by remote-sensing satellite measure- 
ments was first suggested by Kaplan (1959) , and has been the 
object of numerous research studies. 

The Nimbus 3 and Nimbus 4 satellites were equipped with 
SIRS-A and SIRS-B respectively. These first generation SIRS 
instruments, SIRS-A with a vertical scan, and SIRS-B with an 
additional side scan, measured monochromatic radiances in seven 
narrow regions or channels of the 15 ym band of carbon dioxide, 
and one in the window region of water vapor centered at 
899 cm "*■ . Figure 1 (from Gelman et al. 19 72) shows the 
weighting functions dx^/d(ln P) , after Smith et al. (1972), 
of the radiometric emitting channels of the SIRS-B instruments. 
A peak in dx^/dx relative to the In P axis depicts the layer 
in the vertical contributing the maximum radiance to the 
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satellite channel sensor. These radiances could then be 
related to temperature profiles through the Radiative Transfer 
Equation (RTE) 



Here, 1^ is the spectral radiance in channel i (i=l ,2 , . . . , 8) , 
B. [T(P)] is the well-known Planck function for wave number v. 

X -L 

and at temperature T(P), x^(P) is the atmospheric transmittance 
averaged over wave number i (where i signifies the channel 
number centered at v^) , and P is the pressure-level for which 
t^(P) has been computed. 
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FIGURE 1. Derivative of transmittance (x) with respect 



to x(P) = P 2/>7 for SIRS-B 
(after Gelman et al. 1972) 
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Retrieval of temperature profiles from the SIRS-B radi- 
ometer was continued in systematic fashion by the National 
Environmental Satellite Service (NESS) until the radiometer 
output was considered non-usable on 8 April 1971. The number 
of retrievable satellite soundings collected during 1969-71 
was insufficient to make a marked impact on the analysis 
problem over the oceans . 

Anticipating an increased demand for quantitative weather 
observations from orbiting satellites, the National Oceanic 
and Atmospheric Administration (NOAA) in conjunction with the 
National Aeronautics and Space Administration (NASA) , developed 
a second generation set of remote-sensing radiometers, improved 
over the SIRS-type, that would provide twice daily coverage 
of most of the globe and a higher resolution at individual 
scan spots. 

As a result of these developments, the first of the new 
instruments was the Infrared Temperature Profile Radiometer 
(ITPR) as used on Nimbus 5. This radiometer carried four 
carbon dioxide channels, two window channels at 4.3 ym and 
11.5 ym, respectively, and a water vapor channel near 20 ym. 

In 1972, NASA launched the NOAA-2 satellite, which carried 
aloft a new remote- sounding radiometer for which the ITPR 
had been a prototype. The new instrument, the Vertical 
Temperature Profile Radiometer (VTPR) , had been still further 
improved over the ITPR in its angular-array scanning technique. 
The VTPR carried six carbon dioxide channels, a window channel. 
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and a water vapor channel. Observations made at scan spot 
arrays in channels 1,2,... ,8 made it possible to systematically 
retrieve clear-column radiances at a higher resolution than 
before. The new instrument was designed to provide two tem- 
perature profiles per grid point per day on a 400 kilometer 
spacing over the oceans, or a global total of 6000 temperature 
profiles per day. An actual yield of 72%, or about 4000 
profiles per day has been realized [Jastrow and Halem (1973) ] . 
The use of these VTPR clear-column radiances from the NOAA 2 
satellite and the temperature retrieval system described by 
McMillan et al. (1973) was the basis for a thesis by Moran 
(1974) . The basic differences between the Moran study and 
previous work by Smith et al. (1972) was the use, by Moran, 
of standard atmospheric carbon dioxide transmittances with no 
interactive effects on the channel 1,2,..., 6 transmittances 
from the other atmospheric constituents. Moran also computed 
radiances by the use of a 17-layer quadrature scheme for inte- 
gration of Eq. 1 rather than by means of the more usual 100- 
level atmosphere. This was done in order to attempt to reduce 
the trunction error involved with uncertainties in the integra- 
tion of the function [T (P) ] dx^/dp when a relatively crude 
model for x^(P) was utilized. 

The first objective of the present thesis was to improve 
upon the Moran retrieval technique by incorporating recently 
updated state-of-the-art methods [Fritz et al (1972) ] for 
computing atmospheric transmittances, including the combined 
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effects of water vapor, ozone, and carbon dioxide in each 
channel, as well as their vertical path dependencies on 
pressure and temperature. It was also necessary to correct 
these transmittances in each channel for the applicable zenith 
angles of the clear-column spot. For this purpose, a special 
data set for March 1973, including the tuned transmittance 
coefficients for the period, coupled with the first-guess 
temperature and humidity profiles, was acquired from NESS. 

A second and equally important objective was to test the 
final temperature retrievals by comparison with radiosondes 
at island and coastal stations matching the VTPR scan spots 
in both time and spatial location. The retrieved profiles 
were compared at each of the 15 standard pressure levels 
(Table 3) with the matching radiosondes and were analyzed 
both by latitude and by pressure level. 

The retrieval system used here has the feature of inte- 
grating Eq. 1 over 33 atmospheric layers at each clear-column 
radiance scan spot. This system was incorporated in order 
to minimize the truncation error in the quadrature scheme 
applied to the function B^dx^/dp which appears in the RTE 
(Eq. 1). The improved accuracy in the determination of the 
temperature profile using this method is one of the main 
thrusts of this thesis. 
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II. SATELLITE DATA 



The NOAA series of satellites are modifications of the 
earlier Environmental Science Service Administration (ESSA) 
satellites. The major change was the elimination of all 
vidicon camera systems in favor of radiation-sensing instru- 
ments, and in particular the inclusion of the vertical 
temperature profile radiometer (VTPR) . 

A. NOAA 2 DESCRIPTION 

The NOAA 2 satellite orbits the earth at an altitude of 
1464 kilometers once during each 115 minutes. Figure 2 
illustrates the projection of seven orbits onto the earth. 

The solid lines indicate the nortn-to-south or "descending" 
portions of the orbits, which occur over the sunlit hemisphere; 
dashed lines are nighttime "ascending" portions of the same 
orbits. Equatorial crossings occur at 0900 (descending) and 
2100 (ascending) local solar time. 

The VTPR instrument scans perpendicularly relative to the 
subsatellite track from left to right., when viewed downtrack, 
in 23 individual non-overlapping steps covering ±30.3 degrees 
from local nadir. The instrument observes each spot in eight 
spectral intervals or channels through the use of the automatic 
optical filter selector equipment of the VTPR. Six of the 
channels are located in the 15 ym band of carbon dioxide, one 
in the 19 ym region of the rotational band of water vapor, and 
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Figure 2 — Satellite tracks for NOAA 2 VTPR coverage 



one in the 12 ym atmospheric window (see Table 1) . After 
the entire set of 8 spectral measurements has been obtained, 
the scanning mirror rotates to the next spot. 

Table 1. VTPR channel designators and centroid wave numbers 
for filter set four of NOAA 2 VTPR instrument 
number 1 (used in March 1973) . 

CHANNEL 12345678 

WAVE NO. 668.5 677.5 695.7 707.2 724.7 747.5 534.5 835.1 

(cm -1 ) 

B. CLEAR-COLUMN PROCEDURE 

In order to obtain an accurate temperature retrieval, 
the "raw" radiances, that is the radiances measured at the 
VTPR instrument, must be adjusted to eliminate the effects 
of any clouds that may be located in the field of view. 

Smith (1970) and others introduced cloud layer models wherein 
the extent and thickness of cloud layers which would give the 
radiance as sensed at the satellite are mathematically intro- 
duced into the temperature profile calculation. However, the 
results in the neighborhood of a single scan spot tend to 
result in non-unique cloud heights and amounts which causes 
uncertainty in the deduced clear-column radiances. The 
method described below was later developed by Smith (1971) 
in order to eliminate these uncertainties. 

The National Environmental Satellite Service (NESS) 
eliminates cloud contamination through the use of a conversion 
program called CLRAD , applied to a large array of contiguous 
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raw radiance spots. Sets of eight scan lines, containing 
23 scan spots are divided into 3 subarrays or boxes of 8 by 
8, 8 by 7, and 8 by 8 scan spots as shown in Figure 3. Ob- 
servations of raw VTPR radiances are made in each of the eight 
channels and at each scan spot. Figure 4 shows the lateral 
scan geometry along each scan line, where spot 5, spot 12, 
and spot 19 in line 4, are the center of box left, box center 
and box right, respectively. Each box-center is denoted by 
an X in Figure 3. At these locations it is necessary to 
determine the exact zenith angles to the sensor aboard the 
satellite using the following relationship 

SIN ( Z ) = SIN (N ) ^ (2) 

n n z 

Here Z n is the zenith angle at scan spot n (n = 5,12,19) . 

N^ is the nadir angle at the same spots, a is the constant 
6370 km (earth's radius), and H is the altitude of the 
satellite (mean value of 1460 km) . The mean zenith angle at 
spot 12 is 0°, while at spots 5 and 19 it is 23°47' left and 
right of center respectively. 

In addition to the VTPR radiances, measurements in a 
Window-channel Scanning Radiometer (SR), sensing in 10.5 to 
12.5 pm, permits determination of the sea surface temperatures 
at each of the scan spots. McMillin et al. (p 25-33, 1973) 
outline a process of extrapolation of raw-radiance to clear 
column radiance within localized sub-areas about the scan 
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LATITUDE 



ORBIT 




LONGITUDE 



Figure 3 — VTPR scan pattern and data analysis array 
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N = Nadir angle 
Z = Zenith angle 
G = Geocentric angle 
a = Earth's radius (6370 km) 

H = Satellite altitude (mean = 1460 km) 



Figure 4. NOAA-2 VTPR lateral scan geometry at raw 
radiance spots 
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spots. The specifics of the clear radiance conversion method 
involves: (1) the specifications of a reasonably accurate 

value of the first guess temperature profile T(P) to 10 mb, 

(2) an estimation of the water vapor transmittances in the 
window channel derived from the T(P), and (3) the determina- 
tion of the sea surface temperature derived from the SR. The 
value of the clear-column window channel radiance can then be 
computed and compared with the observed window channel radi- 
ance. If the observed "raw" radiance is equal to or greater 
than the computed radiance, the sounding is considered to 
be cloud free. If the computed radiance is greater than the 
"raw" radiance, the sounding is considered to be cloud 
contaminated and the extrapolation process is necessary. 

To facilitate the extraction of the clear-column radiance 
set, "raw" radiance values from adjacent spot pairs (spots 1 
and 2 with their associated I^(v) and I 2 (v) for example) are 
extrapolated to I c i ear ( v )* This can be done if the cloud 
amounts vary linearly in the field of view and share approxi- 
mately the same cloud top level. Figure 5 shows the patterns 
used by NESS for pair comparisons. Each 8 by 8 box, so 
analyzed, will contain 49 scan spots, which can be compared 
with four adjacent spots for a total possible yield of 196 
values of clear-column radiance I c i ear ^ v i^* Here i denotes 
all channels other than the window channel since the window 
channel gives rise only to a possible 49 clear-column esti- 
mates. A proportionately smaller number of both radiance 
sets are found in the 8 by 7 box. If a suitable number of 
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the scan spots are capable of giving, four adjacent dear- 
column values without excessive variances, all resulting 
values are usable in obtaining the weighted clear-column 
radiance in the respective channels. 



STEP BY STEP 
TO LINE 7 




SPOT 



Figure 5. Scanning patterns used by the National Environmental 
Satellite Service (NESS) for pair comparison in 
determining VTPR clear-column radiance. 
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Research on refining the accuracy of determination of 
clear-column radiances is actively being conducted. This 
involves both instrumental accuracy improvements as well as 
the algorithms for the extrapolations of cloud amounts to 
zero. 
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III. TRANSMITTANCE MODELS IN THE VTPR CHANNELS 



A. ATMOSPHERIC TRANSMITTANCE COMPUTATIONS: GENERAL 

CONSIDERATIONS 

Atmospheric transmittances for the six carbon dioxide 
spectral intervals used in the retrieval are fundamental and 
must be computed accurately. The resultant transmittance in 
any given channel and at any pressure level P(J) in the atmos- 
phere may be expressed as the product of the transmittances 
due to the carbon dioxide alone t[C(J)] and the overlapping 
transmittances due to the integrated ozone, OZ(J), and water 
vapor, U(J) , from the top of the atmosphere to the level J 
respectively. Thus 

x(j) = T[C(J) ]*T[U(J) ]*t [OZ(J) ] (3) 

Here C(J) denotes the total absorber mass of the carbon 
dioxide from P(l) = 0.01 mb to P(J) . The symbol OZ(J) de- 
notes the corresponding absorber mass of ozone expressed in 
NTP standard amounts per square centimeter column above level 
J, and U(J) is the precipitable water vapor in the same column. 
As an example of the transmittances produced by the program 
discussed in this thesis, sample outputs are included following 
the program listing. 

Moran (1974) used only carbon dioxide transmittances 
acting in the vertical column (i.e., t [U(J) ]=t [OZ (J) ]= 1, by 
assumption) and further simplified the computation by selecting 
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a mid-latitude, standard-atmosphere • transmittance model for 
carbon dioxide in the six listed channels. Moran found that 
super-adiabatic temperature lapse rates tended to occur 
between the 1000 mb and 900 mb levels. Consequently he found 
it necessary to adjust the values of the transmittance profile 
somewhat arbitrarily in channels 5 and 6, since the trans- 
mittance from the surface in these channels tended to be 
somewhat higher than with the more realistic product 
transmittances of Eq. 3. 

Transmittances in the carbon dioxide channels were based 
upon the laboratory observations of the carbon dioxide infra- 
red spectrum and resultant calculations by Drayson (1966) , 
made on a line-by-line basis assuming a homogeneous path 
(for a fixed temperature) . 

In calculations applied to non-homogeneous vertical paths, 
Smith (1969) , calculates the carbon dioxide transmittances 
corresponding to the U. S. Standard Atmosphere (1962) ; and 
makes improvements on the usual Curtis-Godson modifications 
[Goody (1964)] applied to path transmittances for varying 
pressure, temperature, and absorber masses. In addition, 
modified profiles corresponding to temperature differences 
from the U.S. Standard Atmosphere by ±10, ±20, and ±30°K 
were also computed by Smith, after Drayson' s (1966) laboratory 
calculations. A method for applying temperature corrections 
to the absorber mass in each non-homogeneous sub-layer was 
devised so that the effective temperature-modified carbon 
dioxide transmittances could be adjusted by a. variation of 
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the Curtis-Godson method applied to successive layers, using 
the laboratory data of Drayson. 

The method utilized here for carbon dioxide temperature 
adjustment along the path will be discussed in greater detail 
following that for water vapor. This latter method is some- 
what more general due to the variable mixing ratio of this 
constituent in the vertical as well as the greater complexity 
of the individual line structure in the water vapor infrared 
spectrum. 

B. WATER VAPOR - MIXING RATIO ESTIMATION PROCEDURE 

In order to determine the transmissivity due to water 
vapor in each of the VTPR channels, the precipitable water 
vapor integrated from J = 01 to any arbitrary level must be 
determined. For this purpose an estimation procedure for 
determining the mixing ratio , W(J) , at each level was used 
to represent the water vapor profile. The parameter W(J) 
was used because it was found to be determined from T(P) 
alone with reasonable skill, using the method of Weinreb 
and Crosby (1973) ; and also because it enters directly into 
the calculation of the water vapor transmittances in the 
eight VTPR channels. 

The calculation of the mixing ratio W(J) may be derived 
from 



W(J) = SAT (T, J) *H (J) 



(4) 
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H(J) is the relative humidity and is an unknown function of 
pressure. The use of the water vapor channel (v = 535 cm 
was not attempted for retrieval in this thesis. Furthermore, 
SAT(T,J), the saturated mixing ratio is a function of 
temperature T(J) and pressure P(J) only, and is calculated 
by the GOFF-GRATCH formula [List 1963] 

log e w = -7.90298 (T s /T - 1) + 5.02808 log(T g /T) (5) 

- 1.3816xl0~ 7 [10 11,344 (1 " T/Ts) - 1] 

+ 8.1328xl0 -3 [ 10~ 3 ' 49149 (T s/ T " - 1] + log e 

y ws 

Here e w is the saturation vapor pressure over a plane water 
surface , 

T is the temperature in °K, 

T is the steam point temperature 373.16 °K, and 

e is the saturation vapor pressure of water at the steam 

point temperature . 

The value of SAT(T,J) is then computed using the relationship 

e 

SAT (T , J) = 0.622 (^- rn — ) (6) 

P(J) - e w 

Since no equation directly relates the relative humidity 
to the temperature profile, the mixing ratio was estimated 
from the saturated mixing ratios at a specified set of eleven 
levels (Table 2) by a least-squares formulation based solely 
upon these eleven saturated mixing ratio predictors after 
Weinreb and Crosby (1973) . The method formulates the mixing 
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ratio W(J) at the 40 levels J = 61 , 62 , . . . , 100 and is given 
in the form 



11 

W(J) = W(J) + Z CREG(I,J) (SAT(I) -T(I)] (7) 

1=1 



where W(J) is the sample mean mixing ratio at level J over 

1100 sample soundings that were used to determine 
the regression coefficients, 

CREG(I,J) is a set of predetermined regression coefficients 
calculated for each J-level from the 1100 sample 
soundings , 

SAT (I) is the saturation mixing ratio at the level I, 
where I = 1,2,..., 11 is a subset of the 40 
J-levels that is found to give the best results 
(see Table 2) , and 

T ( I) is the mean saturation mixing ratio at level I 

from the 1100 test soundings. 



Table 2. Eleven levels used for predictor levels in 
water vapor retrieval 



P(J) mb 



1 . 


63 


209.94 


2. 


66. 


245.37 


3. 


70 


299.01 


4. 


77 


412.26 


5. 


82 


509.93 


6. 


90 


699.03 


7. 


92 


753.17 


8. 


94 


810.25 


9. 


96 


870.35 


10. 


98 


933.57 


11. 


100 


1000.00 



Thus the method uses an eleven level saturated mixing 
ratio profile, SAT(I), and converts it into a 40-level mixing 
ratio profile, W(J) , J = 61 , 62 , . . . , 100 . The total 100-level 
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profile is then converted from the 40-level by using the cubic 
power-law profile due to Smith (1966) 

W(J) = W (61) i^l§Y)] 3 (8) 

J<61 

The mixing ratio is then used in determining the water vapor 
transmissivity . 

C. ADJUSTMENT OF WATER VAPOR MASS IN A NON-HOMOGENEOUS 
ATMOSPHERE 

Assuming that the transmittance by water vapor, x[U(J)], 
in a spectral interval of the VTPR is a known function of the 
quantity of the precipitable water vapor U, temperature T, 
and pressure P, an accurate approximation to t[U(J)] in the 
form of a least squares fit polynomial can be calculated. 
Weinreb and Neuendorfer (1973) proposed such a best-fit model 
which in effect becomes an analytic function for calculating 
homogeneous path transmittances . In the Weinreb-Neuendorfer 
method, an inhomogeneous atmosphere is treated as a sequence 
of homogeneous layers, each having a constant temperature, 
pressure, and absorber mass. The transmittance model provides 
for a homogeneous path as a known function of T, P, and U, 
where U is the effective water vapor mass. If t(J) is the 
transmittance in the spectral interval of the instrument at 
level J, then 

t(J) = T [P(J) ,T(J) ,U(J) ] . (9) 
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The procedure involved in defining x [ (J) ] at each level 
as depicted in Fig. 6, uses the functional expression, Eq . 15, 
for layer 1 without modifying the precipitable water vapor 
mass, U, from its initially specified value based upon the 
homogeneous values P(l), T(l), and U(l). This gives rise to 
the following relationship for layer 1 

T ( 1 ) = [P (1) ,T(1) ,U (1 ) ] . (10) 

U(l) = U(l) 



LAYER LEVEL 

1 



3 

n -1 
n 

m (SURFACE) 

Figure 6. Multi-level atmosphere showing P(J), T(J), and 
U (J) structures. 

In proceeding to level 2 , the transmittance through the second 
layer is affected by the new homogeneous characteristics P(2), 
T(2) , and its water vapor increment Al^ , and is no longer along 
the t [P ( 1) ,T (1) ,U (1) ] transmittance curve defined by the 
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conditions in layer 1 (points A^-B^ "^ n 7) • It is 

necessary to determine the adjusted value of water vapor mass, 
denoted V in Fig. 7, with which the radiant transmittance 
enters the new homogeneous layer 2. This is accomplished by 
backing up on curve 2 to U = V by solving 

x (1) = x [P (2 ) ,T (2 ) ,V] (11) 

for V. The solution for V is accomplished by the Newton- 
Raphson iteration. 

This corrected starting water vapor mass, V, with which 
to enter layer 2, has been deduced using curve 2 in Fig. 7. 

For the purpose of computing the transmittance at level 2, the 
new total water vapor mass is then taken as 



W = V + AU 2 (12) 

where AU 2 is the known increase in precipitable water vapor in 
layer 2. This moves the transmittance value x (P (2) ,T (2) ,W] 
down to point B 2 on curve 2 where the value of x(2) is then 
taken as 

x (2) = x[P(2) ,T(2) ,U(2) =W] (13) 

In proceeding downward to level 3 through the next homogeneous 
layer of characteristics P(3), T(3), AU^/ the corrected amount 
of water vapor at level 2 is revised from U(2) to V which 
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T 




I 




V u, w 



Figure 7. Transmittance curves for atmospheric layers 1 and 2 
illustrating modification to effective U(J) in a 
non-homogeneous atmosphere. 
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becomes the level 2 starting point of water vapor mass. 
Hence the new function 



T(3) = t[P(3) ,T(3) , U ( 3) ] 
U ( 3) = V + AU 3 



( 14 ) 



is solved for the new starting mass at the top of level 3. 

The procedure continues as a series of Newton-Raphson itera- 
tions to correct the effective water vapor masses for the 
radiant transmittance upon succeeding layers. This procedure 
is followed through each level J = 2, 3,..., 99 bounding the 
top of the 98 layers below J = 2. 

The procedure is shown by Weinreb and Neuendorfer to be 
exact for a Beer's-Law monochronatic transmittance function. 

D. WATER VAPOR TRANSMITTANCE MODEL 

The water vapor transmittance calculations use the method 
of Weinreb and Neuendorfer (1973) applied to a modification 
of the polynomial representation of water-vapor transmission 
of Smith (1969) . The basic equation for the line transmit- 
tance of water vapor implies that 

14 

In {—In t [U ( J) ] } = Z C T X (15) 

1=1 
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Here 



X 1 


— 


1 


X 8 


— 


X 4 X 7 


X 2 


= 


UT 

(0.1) ln^ 


X 9 


= 


X 3 X 4 


X 3 


= 


In (P/1000) 


X 10 


= 


X 2 X 7 


X 4 


= 


In (T/273) 


X 11 


- 


X 4 X 7 


X 5 


= 


X 2 X 3 


X 12 


= 


(x 4 ) 


X 6 


= 


X 2 X 4 


X 13 


= 


X 3 X 6 


X 7 


= 


(x 2 ) 2 


X 14 


= 


X 3 X 7 



and Cj is a set of coefficients that were predetermined by 
fitting Eq. 15 by the least-squares method using line-by-line 
calculations of the transmittances for 130 homogeneous (con- 
stant pressure and temperature) paths. From Eq. 15 it follows 
that 



t[U(J) ] = exp [-exp ( Z C T X ) ] 

1=1 1 1 

The values of through X^ are directly calculated only 
if U is known. However, X 2 must be calculated indirectly 
since U, the water vapor absorber mass, is essentially a sum 
of scaled values of AU between successive J-levels. Newton's 
method was used to solve for the new value of X 2 / and subse- 
quently other X^ = X ^2 • using the first iteration of 
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+ AU = as the first guess to arrive at the effective 

( 2 ) 

U 2 at level 2, etc. on to the subsequent J-levels. 



Water vapor transmittances in the infrared region have 
been found to contain a spectral feature not previously ex- 
plainable by the classical line-spectral approaches. Bignell 
(1970) found a continuum absorption from the edge of the 20 
ym rotational band which overlaps the lines in the water-vapor 
and carbon-dioxide absorption regions. The continuum absorp- 
tion coefficient was found to be due to the presence of the 
water vapor dimer (t^O^. The dimer absorption coefficient in 

the various spectral regions was found to be proportional to 
-5 

[T(J) ] (see Eqs. 17 below). Thus, it is necessary to add 
the dimer correction to Eq. 15 to fully account for the com- 
bined transmittances of water vapor in all channels. The 
dimer contribution at level J to the ln[-lnx(J)] of Eq. 15 
due to water vapor is given by the adjusted dimer mass profile 
in terms of 



The actual transmittance due to the effective dimer mass 
in each channel is then given by 



DIMER (J) = DIMER ( J-l ) + [ 



A ( J-l) +A ( J) 
2 



] [P (J) -P (J-l) ] 



(17) 




In [-In x (DIMER)] = B (I) *DIMER ( J) 



(18) 
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where B(I), I = 1,2,. ..,8 is a set of 8 constant absorption 
coefficients, the values of which have been given for each 
channel after Bignell (1970). The combined water-vapor 
transmittance in each channel then becomes 

x [U (J) ] = exp{ -exp [DIMER (J) *B(I) ] ) t(J) (19) 

where x (J) is the value computed using Eq. 16. 

The zenith angle correction for water vapor transmittance 
is included in DIMER(J) and U(J) of Eqs . 9 and 17 by making 
the following path-corrections for mass: 

DIMER ( J) = DIMER (J) *PATH, 

U(J) = U ( J) *PATH , (20) 

PATH = sec Z . 

Here sec Z is the secant of the zenith angle from the point 
representing the center of the box to the radiometer (see 
Figure 4) . 

E. CARBON DIOXIDE TRANSMITTANCE MODEL 

A polynomial expression for the carbon dioxide trans- 
mittance x [C(J) ] in terms of logarithmic functions similar 
to Eq. 15 has been shown to be feasible by Smith (1969) . 
However, the absorber mass of carbon dioxide is proportional 
to pressure, and little dependence upon higher powers of the 
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In (P/1000) term has been found necessary to fit line-by-line 
computations of transmissivities with increasing pressure in 
non-homogeneous atmospheres. 

However, the carbon dioxide polynomial expansion, similar 
to Eq. 15 for the water vapor model, has been found to be 
subject to path correction primarily for the temperature varia- 
tion along the vertical. The method of Weinreb and Neuendorf er 
(1973) was shown by Fleming (1974) to be converted into the 
following form 



T [C ( J) ] = t[C(J- 1)] [a(J)+3(J)AT(J)+y(J)AT 2 (J) ] 
+ e i {r[C(J-l) ] — r* ( J— 1) } 



( 21 ) 



where r[C(J)] is the transmittance due to carbon dioxide at 
level J: 

AT(J) is the temperature difference between the U. S. 
Standard Atmosphere (1962) and the average 
temperature over the layer bounded by the 
pressure levels J and J-l (J = 2 , 3 , . . . / 100) , 

a,B,y are a set of 3 by 100 transmittance factors that 
are tuned specifically to the filters of the 
instrument in use for the data period, 

0^ is a set of six adjustment factors applicable to 

systematic corrections for temperature varia- 
tions from the initial standard, and 

T * ( J) is a correction term approximated after the 
Newton-Raphson method from 



T* ( J) 



A ( J) +B ( J) AT ( J) 
AT (J-l) 



D ( J) [AT ( J) ] 
[AT (J-l) ] 2 



( 22 ) 
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where for the initial level J = 1, 



A, B, D are defined as 



A ( 1) = a (1) 

B (1) = B (1) AT (1) (23) 

2 

D(1) = [ ~~ 1 ~ aflT ' + y(1) ] [ AT ( 1) ] 

and for subsequent layers, J > 1, the values of A, B, and D 
are updated iteratively by 



A (1) (J) = A (2) (J-l) 



(2) 



„(1) /TN _ B v (J-l) AT(J) 
B (J) AT (J-l) 



(24) 



D d) (J) = d ( 2) (J-l) [AT ( J) ] 2 

[At(j-i) ] 2 



(2) .... _ D (1) (J) [AT (J) ] 2 a ( J) . * (1) 



D (J) = 



[AT (J-l) ] 



+ A ( J) Y ( J) [AT (J) ] 



B (2) (J) =B (1) (J)a(J) + A (1) (J) B(J) AT(J) 



A (2) (J) = A (1) ( J) a ( J) 



all of which must be computed in the given order. 

F. OZONE TRANSMITTANCES 

The ozone transmittances x[OZ(J)] are minor contributions 
appearing in the form of Eq. 3 for each carbon dioxide channel. 
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They are constant in each channel at a given level, varying, 
however from level to level. Thus the effect of t [OZ (J) ] 
is expressible in the form of a 6 by 100 array, one member 
of the array for each carbon dioxide channel at each J-level. 
The channel values of t [OZ ( J) ] used here are listed in Appendix 
B., as extracted from the NESS CLRAD program. 

The T [OZ (J) ] values in Appendix B have not been path 
corrected and are based upon globally averaged NTP absorber 
paths of ozone in the carbon dioxide channels. 

G. TRANSMISSIVITY ADJUSTMENT FOR ZENITH ANGLE 

When the transmittance in the spectral interval of the 
VTPR has been calculated, an adjustment must be made to 
account for the path differences that occur with various 
zenith angles. The VTPR is designed to measure clear-column 
radiances from points which subtend a zenith angle of 23°47' 
left, 0°, or 23°47' right of the satellite track. These 
angles correspond to the center of the three boxes shown in 
Figure 3. The adjustments to the calculated transmittances 
must be made for these basic angles and any small deviations 
from them. 

The principle of the correction stems from the fact that 
if the zenity angle is other than zero degrees, the atmospheric 
absorbers will have more effect due to the geometry involved 
(Figure 4) with the longer path, and, therefore, the trans- 
missivity in the layer will be decreased level for level. 
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The correction for both the carbon dioxide and ozone 
transmittance at level J is taken from a section of the NESS 
CLRAD program [McMillin (1974 ) ] and is in the form of 

t(J) = t(J) A (25) 

A = (PATH) B (J) (26) 

where PATH is the secant of the zenith angle, and B(J) is a 
polynomial function for carbon dioxide alone. For the ozone 
path corrected transmittances B(J) = 1 for all atmospheric 
levels. B(J) is a function of ln[P(J)] with terms in ln[T(J)] 
up to quadratic (and a product of these logarithmic functions) , 
with a set of five predetermined constants for each channel. 
These constants have been tuned for the particular VTPR 
instrument in use during the period. 
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IV. RETRIEVAL OF TEMPERATURE PROFILES 



A. INPUT DATA USED FOR TEMPERATURE RETRIEVAL 
1 . First Guess Temperatures 

In addition to the climatological soundings, T(P), 
from the U. S. Standard Atmosphere Supplements (1966) (see 
Appendix B) , the retrieval program requires a first guess 
temperature profile in the box center for each VTPR dear- 
column radiance set. The first guess temperatures are of 
two basic types depending upon the latitude of the sounding. 

For tropical latitudes, from the equator to 18 degrees North, 
radiosonde data is averaged over all latitudes and longitudes 
to provide a single climatological first guess profile. For 
the other northern hemisphere retrievals, the first guess pro- 
file comes from an application of the National Meteorological 
Center (NMC) prognostic model for a period of 12 to, at most, 

18 hours from the preceding complete analysis. The profile 
is extrapolated in time and space to the VTPR clear-column 
sounding location. 

The first guess temperatures, as read into the program, 
give temperatures at the 15 mandatory levels shown by an aster- 
isk (*) in Table 3, from 100 mb to 10 mb. The U. S. Standard 
Atmosphere Supplement climatology value is adjusted for latitude 
and month and then fitted to the first guess profile at the 10 
mb level. The adjustment, defined as first guess (10 mb) minus 
standard (10 mb) , is then added to all values of the standard 
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1 

23 

• 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 



56 level atmosphere with merged climatology and 
15 level (*) first guess field for sample sounding 



PRESSURE (mb) 

0.01 

» 

1.00 

• 

10 . 00 * 

15.00 

20 . 00 * 

25.00 

30.00* 

40.00 

50.00* 

60.00 

70.00* 

85.00 

90.00 
100 . 00 * 
150.00* 
200 . 00 * 
250.00* 
300.00* 
400.00* 
500.00* 

600.00 
700.00* 
850.00* 

900.00 

1000 . 00 * 



TEMPERATURE (°C) 
-87.130 

-04.311 



-42.300 

-45.300 

-49.925 

-53.175 

-56.800 

-53.900 

-63.000 

-65.600 

- 68.200 

-69.250 

-69.600 

-70.300 

-63.400 

-58.200 

-48.400 

-38.900 

-23.600 

-11.800 

-03.500 

4.800 

13.600 

16.267 

21.600 
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atmosphere, above 10 mb, and the profile is expanded to the 
56 levels shown in Table 3. 

The 56 level temperature profile is then expanded to 
the 100 J-levels, where the pressure at level J is given by 

7/2 

P(J) = 0 .01 [1+ (J-l) (0 .26087836) ] (27) 

where J = 1,2,..., 100 and defines the pressure from 0.01 mb 
to 1000 mb. 

2 . Other Input Data 

The remaining input data required, having given the 
VTPR first guess sounding, includes the zenith angle, sea 
surface temperature, and clear-column radiances. This data, 
along with all of the verifying radiosondes and first guess 
temperature profiles to 10 mb was kindly provided by NESS in 
the form of an archival tape. 

B. MATHEMATICAL DEVELOPMENT OF THE TEMPERATURE RETRIEVAL 

Given a non-scattering atmosphere in thermodynamic equili- 
brium, the clear-column spectral radiance at the top of the 
atmosphere in each channel is related to the vertical tempera- 
ture profile and the absorber mass structure by the RTE , 

Eq. 1 . 

The right hand side of Eq. 1 is composed of an earth's 
surface radiation term, and an atmospheric contribution which 
is integrated from the surface of the earth (1000 mb) to the 
top of the atmosphere (0.01 mb). 
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When the 100 J-levels developed in Eq. 27 are used, Eq . 1 
can be rewritten 



J=1 



dx 



< (J) 



I. = B. [T(100) ] Tl . (100) + / (B . [T ( J) ] i w \ JT 

1 1 1 j=ioo 1 S3 — )dJ 

(28) 

The Planck radiance function used in Eq. 28 at level J and 
in channel i is 



B ± [T(J) ] 



(Cl)*v i ' 



C2v . 

exp[ TTJT ] “ 1 



(29) 



where is the wave number in channel i, 

-5 2 -1 -1 

Cl is the constant 1.9061x10 erg cm sec ster , and 
C2 is the constant 1.43868 cm °K. 

In order to reduce the effect of truncation error in the 
integration function of the RTE, and increase overall stability, 
the 100-level atmosphere was reduced to a 33-layer atmosphere 
as shown in Fig. 8. In integrating the RTE over the 100 level 
atmosphere, Simpson's rule was applied over each sublayer K 
centered at the fractional J levels of the K-layered atmosphere 
(Fig. 8) . The main virtue of the Simpson quadrature scheme 
using the weighting factors of 1, 4, 1 over each layer (e.g., 

J = 01, 02.5, 04) , lies in the small discretization error, 

5 

(1.5) /90 times the fourth derivative with respect to J of the 
integrand function of Eq. 28, in substitution of the finite 
summation for the Riemann integration. The determination of T 
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LAYER 

(K) 



1 



2 



3 



3 3 



Figure 8 . 



o.oi MB 



LEVEL 

(J) 

1 



2.5 




4 



55 



7 

8.5 




t 0 



t 



9 7 
9 8.5 



1000 MB 



1 0 0 



Thirty three K-layer atmosphere model used for 
the integration of the Radiative Transfer Equation. 
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at levels J = 02.5 etc., is accomplished by interpolation of 

2/7 

the T(P) using the P pressure levels, defined by Eq. 27. 
Therefore, Eq. 28 can be rewritten 

33 

I. = B. [T(100)]T.(100) + Z B. (K)At (K) , (30) 

K=1 

where 

B ± (K) = |{B i [T(J K -1.5) + 4B i [T(K)] + B i [T (J R +1 . 5) ] } 

(31) 

is the layer-mean Planck function for layer K employing 
Simpson's rule. Moreover At^(k) no w spans three J levels 
(Fig. 8) , 



At- (K) = t. (J„-1.5) - t. (J„+1.5) 

1 1 1 x\ 



(32) 



where J K is the value of J at the center of layer K. 

Using the inversion method of Smith (1970) as modified by 
Martin (1974) , the retrieval of mean layer temperatures is 
determined from clear-column radiance values as follows: 

1. Determine the first guess radiance value in 

channel i, , using Eq. 30. 

2. Form the residual by comparing the clear-column 
radiance with the computed value from Eq. 30, using 



(D = T T (l) 

i i-CLEAR i 



(33) 
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3. If the residual is not close to zero the first 

guess temperature profile (P ) must be improved to 

(2) K 

T (Pj ) in accordance with 



[T (2) (Pj 



)] = 



K 



B. [T (1) (P_ ) ] + R. (1) 
1 J K 1 



(34) 



This tends to move the temperature profile closer to the final 
desired convergent profile. If n steps are required to 
converge the residual close to zero, the final step is 



B [T (n+l) (p } j = B [T (n) (p 
1 J K 1 J K 



) ] 



+ Rf n) 

l 



(35) 



where finally 



r(h) 

— < 10" 4 



I i -CLEAR 



(36) 



When a convergent Planck function is determined for each 

channel, i, and at each level, J, the temperature at that 

level can be. determined by solving Eq. 29 for T(J). The 

value of the Planck function in this case must be a weighted- 

mean Planck using the six values of B^(J), i = 1,2,.. .,6, at 

each J-level. Figure 9 shows the values of transmittance in 

2/7 

each of the six channels with respect to P , that is, with 
respect to the linear J-scale. Figure 10, the first derivative 
with respect to J of the transmittances in Fig. 9, is the 
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C0 2 TRANSMITTANCES 




Figure 9. Carbon dioxide transmittances channels 1 through 6, 
after McMillin et al. (1973. 

weighting function used in determining this weighted mean 
Planck. The usual procedure as used by Smith et al. (1972), 
involves the use of a constant normalized wave number 
(v = 700 cm in Eq. 37, whereas in this thesis, a normalized 
wave number, decreasing with decreasing values of pressure was 
found to give a more accurate result, similar to results of 
Moran (1974) . By using this varying wave number (Table 4) , 
advantage can be taken of the level of maximum contribution of 
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Table 4. Normalized wave numbers and corresponding pressure 
levels used in inverting the Planck function for 
temperature retrieval. 



Pressure P(J) 
mb 



1000 - 351 
350 - 96 
95 - 0.01 



Wave Number (v) 
cm"-*- 



720 

710 

690 



C0 2 WEIGHTING FUNCTIONS 




Figure 10. Weighting functions used in determining a 
weighted mean Planck function. 



51 



that channel to the derived equivalent blackbody temperature. 



The weighted mean Planck function then becomes 

(n) 



B (n) tv,T(J) ] = 



6 

Z 

i=l 



B 



(J) 



Ax i (J) 



6 

Z 

i=l 



Ax i (J) 



(37) 



By using this weighted Planck, the 100 mean temperatures, 
T(J), can be determined by the inversion of Eq. 29 



T ( J) = 



In [- 



C2 v 

ClV 3 + B (n) (v,J) 



B (n) (V,J) 



(38) 



In Eqs. 37 and 38, which are accomplished after the iteration 
procedure in each channel is complete, the values of At^(J) 
are centered differences taken over A J = 2 , but centered on 
each individual J level. A forward At^ difference is used to 
initiate the process at level J = 1 and a backward difference 
is used to terminate the weighting process at J = 100. In 
this solution procedure, Eqs. 37, 38, the algorithm used is 
essentially after Smith et al . (1972) , except for the use of 

the variable v. 

C. STEPWISE TEMPERATURE RETRIEVAL 

The procedure for retrieving the layer-mean temperatures 
for layer J = 1,2,..., 100, from the dear-column radiances 
can be summarized as follows (program names refer to the 
listing at the end of this thesis) . 
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1. Input the necessary data. 

a. sounding latitude. 

b. zenith angle. 

c. 15-level first guess temperatures. 

d. sea surface temperature. 

e. clear-column radiances in channels 1 to 6 . 

2. Expand the 15-level first guess temperature 
profile to 56 levels and merge it with climatology at the 
10 mb level using subroutine UPPER. 

3. Expand the 56-level temperature profile to 100 
levels using Eq. 27 and subroutine LEVELL. 

4. Calculate the atmospheric transmittance due to 
carbon dioxide in channels 1 through 6 at each level J using 
subroutine RDTEMP and adjust for the zenith angle by using 
subroutine MARTAU. 

5. Calculate the water vapor content of the atmosphere 
in the form of a mixing ratio profile, one value at each level 
J, using function GOFF and Eqs. 7 and 8. 

6. Use the mixing ratios to calculate the path 

corrected atmospheric transmittance due to water vapor in 
channels 1 through 6 at each level J using subroutine TRANW. 
(Note: steps 5 and 6 were performed on the first iteration 

only, due to the small change in water vapor transmissivity 
from the first guess to the final temperature profile. 



53 



7. Calculate the total atmospheric transmissivity 
in each channel and at each level through the product of 
Eq. 3 which uses path corrected values of transmittance due 
to carbon dioxide, water vapor, and ozone. 

8. For the first iteration, using the 100-level 
first guess temperatures, the Planck functions must be 
calculated in each of the 6 channels. 

a. Calculate the 100 J-level Planck functions 
using Eq. 29. 

b. Calculate the 33 K-level Planck functions by 
fitting a fourth order Lagrangian polynomial to the 100 J-level 
Planck-prof ile [after Conte and deBoor (1965) ] . 

9. Calculate the Ax ( J) for each channel in the 100 
J-level and Ax(K) for the 33 K-level atmospheres. The Ax (K) 
will be used in the discretized form (Eq. 30) of the RTE , in 
order to reduce the truncation error associated with the 
finite summation, substituted for the Riemann integration of 
the integrand function in Eq. 28. This is important while 
dealing with the x(J) values which occur before the temperature 
profile and its associated transmissivities have converged. 

The converged Ax(J) will be used with Eq. 37 in producing the 
new 100-level temperature profile. 

10. Compute the layer-mean Planck values by Eq. 31. 

11. Compute the radiance in each channel using Eq. 30 
and determine the residual difference between the clear-column 
radiance and the computed radiance, Eq. 33. 
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12. If the residual value is not close to zero, 
another iteration must be performed to improve the accuracy 
of the retrieved temperature profile. 

a. Calculate the new 100-level and 33-level 
Planck functions by adding the residual to each previous 
value, Eq. 34. 

b. Calculate the new 100-level temperature 
profile using Eqs. 37 and 38, and go back to step 4. Continue 
recycling until convergence is achieved (see Eq. 36) . 

13. If the residual is small enough (Eq. 36 is true), 
go directly to the calculation of the new 100-level temperature 
profile, Eqs. 37 and 38. 

D. PROGRAM OUTPUT 

The program is designed to produce a mean layer temperature 
profile at 100 levels from 1000 mb to 0.01 mb. For the purpose 
of testing the accuracy of these retrieved temperatures, 
additional inputs were needed. 

First a set of radiosonde soundings, as close as possible 
to each retrieval site in latitude, longitude, and real time 
was read in. While these soundings are not considered to be 
the exact state of the atmosphere at the retrieval site, they 
are as close as possible and therefore were used to determine 
the error in the retrieval program at each level J by using 

ERROR ( J) = T(J) - AT ( J) (39) 
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where T(J) is the retrieved temperature, and 

AT ( J) is the radiosonde sounding temperature. 

Secondly, a set of retrieved temperatures supplied from 
NESS for each case investigated was read in and compared with 
the radiosonde temperatures . 

The program then output the following for each of the 
233 retrievals attempted (see sample output following the 
program listing) : 

1. Station number (1 , 2 , . . . , 233) 

2. Sounding zenith angle in degrees. 

3. Station latitude. 

4. The adjustment at 10 mb necessary to fit the first 

guess temperature profile to the climatological 

temperatures . 

5. The number of iterations necessary for convergence. 

6. Table of results, listing for each of the 15 levels: 

a. Pressure 

b. Retrieved temperature 

c. Radiosonde sounding temperature 

d. Error (retrieved -- sounding) 

e. NESS retrieved temperature 

f. Error (NESS retrieved - sounding) 

In order to easily compare the accuracy of the retrieval 
with the radiosonde data and also to compare the error of 
this method with that of NESS, histograms were drawn and 
statistical data were obtained. The errors were compared 
level for level on an overall basis (globally averaged) and 
also on a latitude-band basis for the latitude ranges shown 
in Table 5. 
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Table 5. Latitude bands used for error analysis of the 
VTPR retrieved temperatures. 



LATITUDE BAND 
1 
2 

3 

4 



BAND WIDTH 
0 £ LAT < 20N 
20N £ LAT < 4 ON 
40N <_ LAT < 60N 
60N < LAT 
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V. RESULTS OF THE VTPR TEMPERATURE RETRIEVAL 



In general, the errors (retrieved temperature - radiosonde 
temperature) were less than 5 °K when averaged over the 
latitude-bands shown in Table 5. The errors also compared 
favorably with the results of the NESS VTPR retrieval in all 
15 atmospheric levels and the four latitude bands. In 
calculating the errors, the assumption that the radiosondes 
matched the VTPR retrieval sites exactly in time and space 
was made. The difference, though small (averaging ±0.6 
degrees in latitude, ±0.7 degrees in longitude, and ±1.87 
hours in time) , coupled with the small data base of 233 
retrievals, could have had some effect on the results and 
should be of some consideration. 

A. GLOBAL RESULTS 

The global mean results obtained here are summarized in 
Table 6 along with the averaged results from the NESS 
retrievals. The frequency distributions of the global error 
analysis at each of 15 key atmospheric levels are shown as a 
function of the error classes encountered (See Appendix A-l 
to A-15) . These frequency distributions, in the form of 
histograms, are included as a visual aid in analyzing the 
deviation of the error from the mean. They show that even 
though the range of errors can extend over as much as 20 °K, 
the probability of a relatively large error is much less than 
the probability of the small errors in the vicinity of the mean. 
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Table 6. Retrieved temperature error analysis on a global basis 



RETRIEVAL ness 



PRESSURE (itb) 


MEAN 

ERROR (°K) 


RLS ERROR 


MEAN 

ERROR (°K) 


RMS ERROR 


NO. OF 
SAMPLES 


10 


-0.55 


4.00 


-0.75 


4.51 


59 


20 


0.57 


3.25 


-0.11 


2.72 


112 


30 


-0.06 


3.20 


-0.17 


3.12 


142 


50 


0.37 


3.20 


-0.17 


3.12 


142 


70 


0.75 


4.39 


0.91 


5.12 


179 


100 


-0.45 


3.89 


-2.30 


5.12 


179 


150 


-0.30 


2.58 


-1.51 


2.73 


208 


200 


-0.66 


3.33 


-1.04 


3.29 


215 


250 


0.15 


3.22 


0.41 


3.33 


206 


300 


-0.58 


2.82 


0.24 


2.70 


223 


400 


-0.18 


2 . 79 


-0.09 


2.07 


226 


500 


-0.67 


3.09 


-0.32 


2.13 


227 


700 


-0.99 


3.29 


-0.86 


2.81 


226 


850 


-0.08 


' 3.26 


-0.31 


3.55 


227 


1000 


-0.24 


1.77 


0.26 


3.89 


196 
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B. LATITUDE BAND RESULTS 



Because of space limitations, frequency distributions 
were not included for the latitude band results. Plowever, 
Table 7 furnishes a summary of the latitude band analysis 
of errors including their root mean square errors along with 
the corresponding NESS results. 
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LATITUDE BANDS IN DEGREES NORTH 
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LATITUDE BANDS IN DEGREES NORTH 
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VI. CONCLUSIONS 



The first objective of this thesis was to incorporate the 
state-of-the-art methods [as used by Fritz et al . (1972)] in 

calculating atmospheric transmittances , into the retrieval 
scheme. The method involved here for calculating the effects 
of water vapor, ozone, and carbon dioxide transmittance proved 
to be satisfactory. This, coupled with the inclusion of a 
33 level integration scheme and the use of a variable wave 
number, as opposed to the constant v = 700 cm \ combined to 
give results comparable to those of NESS. 

The temperature profile convergence occurred after an 
average of 5 iterations. While this is an improvement over 
the Moran (1974) model, as would be expected with the more 
sophisticated procedures incorporated, the number of iterations 
required would possibly be reduced further through the use of 
a more accurate first guess temperature profile. Gelman et 
al . (1972) used a regression technique for estimating the 

first guess temperatures above 10 mb, as opposed to merging 
with the standard atmosphere climatology as used here. Since 
this regression technique of Gelman produces a more accurate 
first guess temperature profile, it would reduce the number 
of iterations required for convergence even further. However 
this program was not available locally for inclusion in this 
thesis . 
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The second objective, that of evaluation through compari- 
son with matching radiosondes, showed that the retrieval 
system produced errors that were comparable to NESS retrievals, 
or better, especially in the area of the tropopause. The 
small data base of 233 soundings is not considered large 
enough to evaluate the accuracy of the retrieval system beyond 
that, and more experimentation with a larger sample, and 
inclusion of the Gelman regression technique is required. 

It is concluded that this thesis has demonstrated a 
consistent and accurate means of temperature profile retrieval 
from VTPR clear-column radiance data. 
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APPENDIX A-l 



Histogram of the globally averaged errors showing probability 
distribution [P (E) ] of the occurrence of the errors, the actual 

J£ r0 M fre m? nCle ! ln the err ° r class E (marked along the top) , 
the Mean (M) , and the Quartiles (•) , for sample size 196. 
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APPENDIX A-2 



Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 227. 
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APPENDIX A- 3 

Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 226. 
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APPENDIX A-4 



Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•)/ for sample size 227. 
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APPENDIX A- 5 

Histogram of the globally averaged errors showing probability 
distribution [ P (E ) ] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•) , for sample size 226. 
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APPENDIX A- 6 



Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•) , for sample size 223. 
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APPENDIX A- 7 



Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 206. 
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APPENDIX A- 8 

Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•)/ for sample size 215. 
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APPENDIX A- 9 



Histogram^ of the globally averaged errors showing probability 
distribution [P (E) ] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 208. 
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APPENDIX A-10 



Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•)/ for sample size 179. 
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APPENDIX A-ll 



Histogram of the globally averaged errors showing probability 
distribution [ P (E) ] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (♦) , for sample size 179. 
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APPENDIX A-12 

Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) 
the Mean (M) , and the Quartiles (•)/ for sample size 142. 
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APPENDIX A-l 3 



£n r X f avera 9 ed errors showing probability 

distribution [P(E)] of the occurrence of the errors, the actual 

the°L5n e< ?M? nCle ^ iu the err ° r class E (marked along the top) , 
t e Mean (M) , and the Quartiles (•), for sample size 142. 
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APPENDIX A- 14 

Histogram^ of the globally averaged errors showing probability 
distribution [P(E) ] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 112. 
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. APPENDIX A-15 

Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
tu r °l fre( 3uencies in the error class E (marked along the top) , 
the Mean (M) , and the Quartiles (•), for sample size 59. 
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APPENDIX B 



OZONE TRANSMITTANCES FOR THE SIX VTPR CLE AR-C3LUMN CHANNELS 
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CARBON DIOXIDE TR AN SMI TTAK'CE S FOR HE SIX VTPR CHANNELS 



LEVEL 


CHAN 1 


CHAN 2 


CHAN 3 


CHAN 4 


CHAN 5 


CHAN 6 


1 


0. 991131 


0.999653 


0.999886 


0.999912 


0.999936 


0.999993 


2 


0.988616 


0.999418 


0.999733 


0. 999799 


0.999842 


0.999935 


3 


0. 985314 


0.999092 


0 . 9°95 3 7 


0.999575 


0.999653 


0.959962 


4 


0. 980754 


0.998608 


0.999259 


0.999173 


0.999303 


0.969909 


5 


0. 975299 


0.997973 


0.998884 


0.998641 


0.998848 


0.999315 


6 


0. 969083 


0.997204 


0.998390 


0.998053 


0.998348 


0.999674 


7 


0.96 1684 


0.996247 


0.997750 


0. 997415 


0. 997822 


0.999486 


8 


0.952678 


0.995045 


0.996952 


0.9° 6 73 8 


0.997281 


0.999253 


9 


0. 941640 


0.993540 


0.995971 


0.996022 


0.996731 


0.998991 


10 


0.928429 


0.991679 


0. 9947 89 


0. 995254 


0.9961 56 


0.958697 


11 


0.913059 


0.989418 


0.993404 


0.994433 


0.995537 


0.998400 


12 


0.895530 


0.9867 19 


0.991792 


0.993556 


0.994859 


0. 998101 


13 


0. 875965 


0.983537 


0.989937 


0.992624 


0. 99 41 12 


0.907302 


14 


0. 854614 


0.979809 


0.987313 


0.991627 


0.993282 


0.997503 


15 


0.831837 


0.975464 


0.985379 


0.990533 


0.992342 


0. 9971 93 


16 


0.808120 


0.970429 


0.982596 


0.989324 


0.991280 


0.996866 


17 


0. 7839 78 


0.964630 


0.979476 


0.987985 


0.990083 


0.996519 


18 


0. 7598 06 


0.957937 


0.975819 


0.986496 


0.938739 


0.996145 


19 


0.735895 


0.950424 


0.97] 736 


0. 984847 


0.987244 


0.995739 


20 


0.712363 


0.941861 


0.967135 


0.983030 


0.935595 


0.995301 


21 


0.689153 


0.932223 


0.961984 


0.931041 


0.983792 


0.994323 


22 


0. 666067 


0.921432 


0.956255 


0.978881 


0.981841 


0.994319 


23 


0.642342 


0.909414 


0.949899 


0.975531 


0.979731 


0.993767 


24 


0.619215 


0.896098 


0.942897 


0.973993 


0.977470 


0.993178 


25 


0. 594999 


0.801417 


0.9352 16 


0.971258 


0. 975059 


0.992543 


26 


0. 570067 


0.865308 


0.926777 


0.953284 


0.972 4 73 


0 . 991 863 


27 


0. 544428 


0.847723 


0.91 7538 


0.95505 7 


0.959713 


0.991121 


28 


0. 518133 


0.828618 


0.907477 


0. 951 575 


0.966787 


0.99032L 


29 


0.491282 


0.807966 


0.896576 


0.957833 


0.963705 


0.989465 


30 


0. 463997 


0.785753 


0.884823 


0.953848 


0.960479 


0. 988555 


31 


0.43641 4 


0.761986 


0.872210 


0.949607 


0.957125 


0, 987592 


32 


0. 408677 


0. 73i>687 


0. 858738 


0.945123 


0.953658 


0.936582 


33 


0. 380939 


0.709907 


0. 8444 14 


0.940408 


0.950097 


0.985527 


34 


0. 353356 


0.681707 


0.8292 54 


0.935474 


0. 946459 


0.984435 


35 


0. 326084 


0.652176 


0.813267 


0.930328 


0.942757 


0. 983308 


36 


0. 299286 


0.621431 


0.796393 


0.924929 


0.938977 


0.982136 


37 


0. 273120 


0.589637 


0. 778590 


0.919237 


0.935107 


0.980904 


38 


0.247735 


0.556957 


0. 759887 


0.913260 


0.931157 


0.979613 


39 


0. 22 32 74 


0.523557 


0. 740321 


0. 907008 


0. 927138 


0.973282 


40 


0.199827 


0.489605 


0.720067 


0.900580 


0.923095 


0.976925 
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41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 



0 . 177535 
0 . 156501 
0 . 136779 
0 . 118460 
0 . 101568 
0 . 086144 
0.072217 
0.059791 
0 . 048 S 30 
0.039287 
0.031120 
0 . 024258 
0 . 018534 
0.013970 
0.010292 
0 . 007424 
0.005227 
0 . 003538 
0 . 002388 
0 . CO 1 54 1 
0. 000968 
0. 000591 
0.000341 
0.000180 
0 . 000082 
0. 000020 
0.000000 
0 . 0 
0.0 
0 . 0 
0.0 
0.0 
0 . 0 
0 . 0 
0.0 
0.0 
0 . 0 
0 . 0 
0 . 0 
0.0 



0.455338 
0.420970 
0.386711 
0.352818 
0.319468 
0.286895 
0.255350 
0.225062 
0.196242 
0.169099 
0.143870 
0 . 120747 
0.099830 
0.081190 
0.064897 
0.050927 
0.039151 
0 . 02943 ? 
0.021629 

r» A i C C i c 

0.010838 
0.007353 
0.004810 
0.003047 
0.001836 
0 .001156 
0.000691 
0.000395 
0.000221 
0.000123 
0.000060 
0.000019 
0.000000 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 



0.699027 
0.677231 
0.654787 
0.631611 
0.607929 
0. 583764 
0. 559096 
0.533982 
0.508515 
0.482333 
0 . 456970 
0.430848 
0.404528 
0.378156 
0. 351884 
0.325859 
0 . 300236 
0.2752 13 
0.250973 
0.227612 
0.205211 
0 . 183853 
0. 163613 
0.144515 
0.126586 
0 . 109356 
0.094406 
0. 080334 
0.057655 
0.056346 
0.046403 
0.037766 
0.030335 
0.023984 
0.01 8661 
0 . 01429 ? 
0.010759 
0. 007967 
0.005844 
0.004220 



0.893880 
0. 886900 
0.879668 
0.872086 
0.864268 
0.856180 
0.847746 
0.838040 
0 . 82 9 7 70 
0. 820248 
0.810294 
0.799742 
0. 783502 
0. 775550 
0.763863 
0.750427 
0.736235 
0 . 721 332 
0.705756 
0.609479 
0.672462 
0. 654664 
0.536031 
0.616461 
0.595843 
0. 574072 
0.551202 
0 . 52744 ? 
0 . 532 S 86 
0.477534 
0.451509 
0.425025 
0. 393223 
0.371034 
0.343825 
0.315940 
0. 290388 
0.254438 
0. 239271 
0.215037 



0.918984 
0.914800 
0.91 0554 
0.906201 
0.901789 
0.897299 
0 . 892698 
0.887974 
0.8931 25 
0. 878152 
0.873021 
0.857663 
0.862036 
0.856124 
0 . 849909 
0.843377 
0 . 8365 17 
0.829334 
0.821833 
0.013992 
0. 805783 
0. 797179 
0. 788142 
0.773612 
0.758520 
0.757793 
0.746438 
0.734486 
0.721936 
0. 703741 
0.694903 
0.680453 
0.665387 
0.649593 
0.633165 
0.616200 
0.598623 
0. 583505 
0 . 56 18 78 
0.542755 



0.975523 
0.974366 
0.972573 
0.971304 
0.969402 
0.967754 
0.966337 
0.964244 
0.962376 
0.960436 
0 . 95&431 
0.956221 
0.953867 
0.951329 
0.943596 
0.945657 
0.942506 
0.939146 
0.935582 
0.931796 
0.927765 
0.923464 
0.918861 
0. 913897 
0.903503 
0.902609 
0.896197 
0.889293 
0.881867 
0.873871 
0. 865283 
0.856134 
3.846389 
0.835938 
0.824852 
0.813236 
0.800919 
0.783343 
0.774598 
0. 760583 
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81 


0 . 0 


0.0 


0.002995 


0 . 191390 


0.523149 


0 .745995 


82 


0 . 0 


0.0 


0.002071 


0 . 170036 


0.503112 


0. 730874 


83 


0 . 0 


0.0 


0.001380 


0.149654 


0.482760 


0.715302 


84 


0.0 


0.0 


0.000385 


0 . 130787 


0.452181 


0.699323 


85 


0.0 


0.0 


0.000548 


0.113420 


0.441413 


0.682954 


86 


0.0 


0.0 


0.000334 


0.097547 


0.420495 


0.666182 


87 


0 . 0 


0.0 


0.000208 


0.083170 


0.399466 


0.649015 


88 


0.0 


0.0 


0.000132 


0.070295 


0.378367 


0.631455 


89 


0.0 


0.0 


0.0 


0.05891 5 


0. 357233 


0 . 61*499 


90 


0.0 


0.0 


0.0 


0.048942 


0.335122 


0.595161 


91 


0.0 


0.0 


0.0 


0.040301 


0.315182 


0.576550 


92 


0.0 


0.0 


0.0 


0. 032900 


0. 294559 


0.557781 


93 


0.0 


0.0 


0.0 


0. 026612 


0.274311 


0.538879 


94 


0.0 


0.0 


0.0 


0.0213 19 


0.254492 


0.519862 


95 


0 . 0 


0.0 


0.0 


0.015910 


0.235156 


0.500744 


96 


0.0 


0.0 


0.0 


0.013286 


0 . 216375 


0.481554 


97 


0.0 


0.0 


0.0 


0.010351 


0.198231 


0.462330 


98 


0 . 0 


0.0 


0.0 


0.008008 


0.180793 


0.443100 


99 


0.0 


0.0 


0.0 


0.005171 


0. 1641 26 


0.423879 


1 00 


0.0 


0.0 


0.0 


0.004732 


0.148023 


0 . 404? 73 



Ill 



WATER VAPOR TR ANSM I TT ANCE S FOR THE SIX VTPR CHANNELS 



VEL 


CHAN 1 


CHAN 2 


CHAN 3 


1 


1.000000 


1.000000 


1.000000 


2 


1.000000 


1.000000 


1.000000 


3 


1. OOOOOO 


1.000000 


1.000000 


4 


1. OOOOOO 


1.000000 


l . OOOOOO 


5 


1. OOOOOO 


1.000000 


1.000000 


6 


1 . OOOOOO 


1 .OOOOOO 


1.000000 


7 


1.000000 


1 .OOOOOO 


1.000000 


8 


1. OOOOOO 


1.000000 


1 . OOOOOO 


9 


1. OOOOOO 


1 .OOOOOO 


1.009000 


10 


1.000000 


1 .OOOOOO 


1.000000 


11 


1.000000 


1.000000 


1.000000 


12 


1.000000 


1 .OOOOOO 


1.000000 


13 


1.000000 


1 .OOOOOO 


1.000000 


14 


1.000000 


1 . OOOOOO 


1.000000 


15 


1 . OOOOOO 


1 .OOOOOO 


1.000000 


16 


1. OOOOOO 


1.003000 


1.000000 


17 


1. OOOOOO 


1.000000 


1.000000 


18 


1.000000 


1 .OOOOOO 


1.000300 


19 


1 . OOOOOO 


1 .OOOOOO 


1.000000 


20 


1.000000 


1 .OOOOOO 


1 . OOOOOO 


21 


1. OOOOOO 


1 . OOOOOO 


1.000000 


22 


1. OOOOOO 


1.000000 


1.000000 


23 


1.000000 


1 .OOOOOO 


1.000000 


24 


1. OOOOOO 


1.000000 


1.000000 


25 


1.000000 


1 .OOOOOO 


1 . OOOOOO 


26 


1.000000 


1.000000 


1.000000 


27 


1.000000 


1.000000 


1.000000 


28 


1. OOOOOO 


1 .OOOOOO 


1.000000 


29 


1. OOOOOO 


1.000000 


1.000000 


30 


1.000030 


1 . OOOOOO 


1 . OOOOOO 


31 


1. OOOOOO 


1.000000 


1.000300 


32 


1.000000 


1 .OOOOOO 


1.000000 


33 


1.000000 


1 .OOOOOO 


1.000000 


34 


1.000000 


1.000000 


1.000300 


35 


1.000000 


1.000000 


1.000000 


36 


1.000000 


1 .OOOOOO 


1.000000 


37 


0.999999 


1.000000 


1.000000 


38 


0. 99999Q 


1 .OOOOOO 


1.000000 


39 


0. 999998 


1.000000 


1.000000 


40 


0. 999997 


1.000000 


1.000000 



CHAN 4 



1 . 030000 

1 .033000 

1 . 003000 
1.033300 

1 . 030000 

1 . 030000 

1.000000 

1 . OOOOOO 

1 . 000300 

1 . OOOOOO 

1 . OOOOOO 

1 .003000 

1 . OOOOOO 

1 . 000000 

1.000000 

1 . 030000 

1 . 030000 

1.033000 

1 . 003000 

1.033000 

1 . 000000 

1 . 003000 

1.000000 

1 . 030000 

1 . 033000 

1 . 030000 

1 . OOOOOO 

1.033000 

1.030000 

1 . 000000 

1.030000 

1 . 030000 

1 . 003000 
l . 033000 

1 . 000000 

1 . 033000 

1.000000 

1 . 003000 

1 . OOOOOO 

1 . 033000 



CHAN 5 



1.000000 

1 . 000000 
1.000003 

1.000000 

1.000000 

1 . 000000 

1.000000 

1 . 000000 

1.000000 

1 . 000000 

1. 000000 

1.000000 

1.000000 

1 . 000300 

1.000000 

1 . 003000 

1 . 000000 

1.000000 

1. 000300 

1.000000 

1.000000 

1 . 003000 

1.000000 

1 . 030000 

1.000000 

1.003000 

1.000000 

1.000000 

1.003000 

1 . 000300 

1 . OOOOOO 

1.003000 

1 . OOOOOO 

1.003000 

1 . OOOOOO 

1 . 000000 

1.000000 

1 . OOOOOO 

1.003000 

1 . 003000 



CHAN 6 



1 . 000033 

1.000033 
1 .003303 

1.000033 

1.000003 

1.000303 

1.000003 

1.000033 

1.000033 

1.000033 
1 .003033 
1 .000033 

1.000003 

1.000333 

1.000003 
1.003333 

1 . 000033 

1 . 000033 

1 .000003 

1.000033 

1 . 000303 

1.000303 

1 .000303 

1.000333 

1 . 000003 

1.000333 

1.000303 
1.000000 

1.000003 

1 . 000333 

1.000003 

1.000003 

1 . 000003 

1.000003 
l . 000030 

1 .000033 

1.000003 

1.000333 

1 .000003 
1.003303 
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41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 



0.999996 
0. 999993 
0. 999990 
0.999985 
0.999978 
0.999968 
0. 999955 
0.999938 
0.999915 
0.999884 
0.999845 
0.999794 
0.999728 
0. 999645 
0. 999540 
0. 999409 
0.999246 
0. 999046 
0. 998802 
0. 998508 
0.998154 
0. 997638 
0.996922 
0.996058 
0.995036 
0. 993836 
0.992429 
0.990770 
0.988810 
0 . 986513 
0 . 983706 
0 . 980191 
0.975928 
0. 970894 
0.965046 
0. 958313 
0.950473 
0 . 941307 
0.930816 
0.919023 



1.000000 

1.000000 
0.999999 
0 . 999^99 
0 . 999999 
0.999998 
0.999997 
0.999996 
0.999994 
0.999991 
0.999987 
0.999982 
0.999975 
0.999965 
0.999951 
0.999933 
0.999909 
0.999877 
0.999835 
0.999782 
0.999714 
0.999609 
0.999451 
0.999246 
0.998985 
0.998659 
0.998254 
0.997752 
0.997129 
0.996367 
0.995397 
0.9941 34 
0.992542 
0.990595 
0.988257 
0.985478 
0.982140 
0.97811 1 
0.973349 
0.967821 



0.999999 
0. 999999 
0.999998 
0.999997 
0. 999996 
0.999993 
0.999991 
0.999987 
0.999981 
0.999973 
0.999963 
0.999949 
0.999929 
0.999902 
0.999865 
0.999316 
0.999750 
0.999663 
0.999550 
0.999406 
0.999221 
0. 998936 
0.998507 
0.997946 
0.997228 
0.996322 
0.995187 
0.993765 
0.991987 
0.989793 
0 . 986988 
0.933340 
0.9787 58 
0.9731 75 
0 . 9665 07 
0.958637 
0.949299 
0.938244 
0.925471 
0.91 1009 



0.999999 
0.999999 
0.999998 
0.999997 
0.999995 
0.999993 
0.999990 
0.999985 
0.999979 
0.999970 
0.999959 
0.999943 
0. 999921 
0.999891 
0.999352 
0.999799 
0.999730 
0.999640 
0.999525 
0. 999380 
0.999196 
0.998915 
0,998501 
0.997970 
0.997302 
0.996477 
0.995462 
0.99421 1 
0.992674 
0.990807 
0.938452 
0.965424 
0.931 663 
0.977128 
0.971763 
0.965490 
0.958100 
0.949395 
0.939377 
0.928071 



1 . 000000 

1 .000000 

1 . 000000 

1.000000 

1 . 000000 

1.000000 
0.999999 
0.999999 
0.999998 
0.999997 
0. 999996 
0.999995 
0.999992 
0.999988 
0.999983 
0.999976 
0.999966 
0.999952 
0.999934 
0.999910 
0.999878 
0.999827 
0.999749 
0.999644 
0. 999507 
0.999331 
0.999107 
0.998823 
0. 993464 
0. 998014 
0.99 1 429 
0.996650 
0.995646 
0.994392 
0.992855 
0.990997 
0.983726 
0.985944 
0.982611 
0 . 978699 



1.000000 

1 .000000 

1 .000000 

1.000000 

1.000000 
0.999999 
0.999999 
0.999999 
0.999993 
0.999997 
0.999995 
0.999994 
0.999991 
0.999986 
0.999980 
0.999971 
0.999958 
0.999940 
0 . 99991 6 

r\ o rv / 

KJ • 7;ooc-t 

0 . 9993^1 
0.999772 
0.999662 
0.999511 
0. 999303 
0.999040 
0.993691 
0.998236 
0.997646 
0.996391 
0.995893 
0.994548 
0.992799 
0.990596 
0.987877 
0.984569 
0.980522 
0.975585 
0.969703 
0.962857 



113 



81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 



0.905727 
0 . 890766 
0. 873758 
0. 854578 
0 . 833367 
0.810018 
0. 783846 
0. 754456 
0.721769 
0.685355 
0 . 6 * 4-4398 
0. 596892 
0. 541886 
0.480873 
0. 416116 
0.349820 
0. 283816 
0.219603 
0 . 160078 
0. 110715 



0.961377 

0.953864 

0.944969 

0.934496 

0.922399 

0.908478 

0.892089 

0.872718 

0.850020 

0.823326 

0.791447 

0.751831 

0.702560 

0.643995 

0.577550 

0.504901 

0.427585 

0.346967 

0.266685 

0.195466 



0.894622 
0 . 8761 60 
0.855273 
0. 831906 
0.806297 
0.778395 
0.747524 
0 . 7133 74 
0.676024 
0.635195 
0. 590386 
0. 540043 
0.483718 
0.423278 
0.361101 
0.299315 
0.239564 
0.183059 
0.132022 
0.090253 



0.91 5290 
0.900897 
0.884572 
0 . 856 2 26 
0. 846008 
0.823828 
0.799060 
0. 771343 
0. 74061 1 
0.705457 
0.563154 
0.623868 
0.572615 
0.515535 
0.454401 
0.390912 
0.326423 
0. 261954 
0.199974 
0.145974 



0. 974090 
0.958669 
0.962200 
0.954532 
0.945613 
0.935296 
0.923055 
0.908470 
0.391226 
0.870735 
0.845954 
0.81 4646 
0.774398 
0 . 72649 ? 
0.659975 
0. 606039 
0.535146 
0.457417 
0.375108 
0. 296695 



0.954860 
0.945573 
0.934722 
0.922169 
0.907931 
0.891857 
0.873369 
0.852050 
0.827689 
0.799774 
0.767475 
0.728875 
0.682605 
0.629166 
0.569725 
0. 505494 
0.437412 
0.365077 
0.293774 
0.226899 
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\ 






TOTAL ATMOSPHERIC TR AN SM I TT A NC E ( PATH-CORRECTED) 



LEVEL 


CHAN 1 


CHAN 2 


CHAN 3 


CHAN! 4 


CHAN 5 


CHAN 6 


1 


0.991123 


0.999640 


0.999872 


0.999899 


0.999927 


0.999985 


2 


0.988599 


0.999392 


0.999705 


0. 999767 


0.999816 


0.999970 


3 


0.985280 


0.999045 


0.999485 


0.999515 


0.999602 


0 . 999933 


4 


0. 980696 


0.998526 


0.9591 73 


0.999073 


0.999220 


0.999957 


5 


0. 975207 


0.997838 


0.998750 


0. 999485 


0.998706 


0.999729 


6 


0. 968946 


0.996999 


0.9981 94 


0.997825 


0.998129 


0.999543 


7 


0.961489 


0.995951 


0.997475 


0.997095 


0.997501 


0.999296 


8 


0. 952409 


0.994633 


0.996579 


0.995305 


0.996834 


0.993992 


9 


0. 941280 


0.992984 


0.995477 


0.995449 


0.996125 


0.998631 


10 


0.927962 


0.990948 


0.994151 


0.994515 


0. 995356 


0.998223 


11 


0.91 2467 


0.988479 


0.992592 


0.993497 


0.994508 


0.997787 


12 


0. 894798 


0.985537 


0.990779 


0.992390 


0.993560 


0.997325 


13 


0. 875073 


0.982071 


0.988691 


C. 991189 


0.992497 


0.996836 


14 


0. 853544 


0.978018 


0.986298 


0.989880 


0.991301 


0.996315 


15 


0.830569 


0.973299 


0.983554 


0.938427 


0.989940 


0.995750 


16 


0. 806635 


0.967841 


0.930417 


0.936804 


0. 938397 


0.995130 


17 


0. 782254 


0.961565 


0.976845 


0.984995 


0.986654 


0.994443 


18 


0. 757824 


0.954394 


0.972786 


0.982973 


0.984699 


0.993697 


19 


0.733641 


0.946252 


0.968199 


0.980733 


0.982529 


0. 992373 


20 


0* 70982 7 


n 0-J-7A 79 

0.’926781 


0.963047 


o m o *> t. o 

y-/ m y 1 6. vO 


0.980144 


0.901973 


21 


0.686330 


0.957293 


0.975577 


0.977551 


0.991000 


22 


0.662957 


0.915307 


0.950914 


0.972658 


0.974753 


0.989950 


23 


0.639437 


0.902565 


0.943862 


0.959483 


0.971724 


0. 9S881 0 


24 


0.615498 


0.888472 


0.9361 19 


0.966039 


0.968449 


0.987580 


25 


0. 590951 


0.872951 


0.927666 


0.962301 


0.964909 


0.986250 


26 


0. 565674 


0.855943 


0.918391 


0.958227 


0.961078 


0. 984804 


27 


0.539702 


0.837437 


0.908280 


0. 953816 


0.956982 


0.983241 


28 


0.513110 


0.817422 


0.897295 


0.949086 


0.952663 


0.981568 


29 


0. 486019 


0. 795908 


0.8854 15 


0.944054 


0.948159 


0.979792 


30 


0.458565 


0.772911 


0. 872623 


0. 933743 


0.943515 


0.977917 


31 


0.430898 


0.748469 . 


0.858910 


0.933172 


0.938775 


0.975952 


32 


0. 403163 


0.722624 


0.844280 


0. 927363 


0.933976 


0.973906 


33 


0.375505 


0.695429 


0. 828747 


0.921333 


0.929138 


0.971791 


34 


0. 348059 


0.666944 


0. 812395 


0.915100 


0.924278 


0.969614 


35 


0. 320999 


0.537259 


0.795206 


0.903677 


0.919405 


0.967333 


36 


0. 294448 


0.606487 


0.7771 59 


0.902027 


0.914505 


0.965107 


37 


0.268565 


0.574790 


0. 758231 


0.895116 


0.909567 


0.962764 


38 


0. 243488 


0.542324 


0. 738474 


0. 897957 


0. 904598 


0.960368 


39 


0. 219354 


0.5092 50 


0. 71 7938 


0.880566 


0.899609 


0.957930 


40 


0. 196243 


0.475728 


0.6963 14 


0.873041 


0.894641 


0.955485 
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41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 



0 . 174290 
0 . 153591 
0 . 134196 
0 . 116190 
0. 099594 
0.084448 
0.070776 
0. 058581 
0.047828 
0.038468 
0.030461 
0.023735 
0 . 0181 76 
0.013657 
0. 010057 
0 . 007251 
0. 005102 
0. 003500 
0. 002329 
0. 001502 
0. 000943 
0.000575 
0.000332 
0. 0001 75 
0. 000030 
0.000019 
0 . 000000 
0.0 
0 . 0 
0.0 
0 . 0 
0 . 0 
0.0 
0 . 0 
o.c 
0 . 0 
0.0 
0.0 
0.0 
0 . 0 



0.441985 
0.408228 
0.374654 
0. 341508 
0.308956 
0.277221 
0.246539 
0.217125 
0.139178 
0.162892 
0.138491 
0.116153 
0.095968 
0.078000 
0.062309 
0.048868 
0.037548 
0 . 0282.13 
0.020723 
0.014858 
0.010375 
0.007036 
0.004601 
0.002913 
0.001803 
0.001104 
0.000659 
0.0003 76 
0.000210 
0.0001 17 
0.000057 
0.000019 
0.000000 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 



0.675019 
0. 652599 
0. 629669 
0.6061 61 
0.582299 
0. 5581 13 
0.533589 
0 . 508783 
0.483785 
0.458732 
0.433634 
0.400421 
0.383130 
0.357879 
0.332799 
0.3030 14 
0.283658 
0. 259906 
0 . 2369 23 
0 . 214790 
0.1935 78 
0. 173351 
0.1541 80 
0 . 136078 
0 . 119093 
0 . 103241 
0. 088606 
0.075280 
0.0632 76 
0.052575 
0.0431 69 
0.035000 
0.027979 
0.021993 
0.016993 
0.012907 
0.009621 
0.007040 
0.005094 
0.003621 



0. 855296 
0.857326 
0.849163 
0. 840715 
0.832097 
0.823279 
0.814191 
0.304812 
0.795153 
0.785227 
0. 774961 
0.764193 
0. 752833 
0. 740355 
0.728234 
0.714949 
0. 700992 
0.685396 
0.671196 
0.653358 
0.538837 
0.521560 
0.603452 
0.584433 
0.554389 
0.543216 
0.520958 
0.497802 
0.473828 
0.449036 
0.423501 
0.397390 
0.370863 
0.343904 
0.315897 
0.290197 
0.263819 
0 . 233032 
0.213081 
0 . 189172 



0. 797749 
0.884633 
0.879594 
0.874495 
0.869376 
0. 864220 
0.858994 
0. 853585 
0.848290 
0.842810 
0.837210 
0.831424 
0.325415 
0 . 8191 69 
0. 812670 
0 . 805906 
0.798868 
0.791560 
0. 783989 
0 . 776134 
0. 767968 
0. 759458 
0.750562 
0.741219 
0. 731356 
0.720895 
0.709321 
0.698165 
9.685906 
0.672989 
0.659386 
0.645095 
0.630102 
0.614297 
0.597764 
0. 580587 
0.562664 
0.544031 
0.524729 
0.504789 



0.953010 
0.950512 
0.948003 
0.945461 
0.942924 
0.940389 
0.937840 
0.935272 
0.932694 
0.930114 
0.927513 
0.924345 
0.922074 
0.9191 85 
0.916160 
0.912987 
0.909653 
0.906147 
O.QO?^77 
0.893614 
0.894529 
0.893181 
0.835520 
0.880433 
0.875310 
0.869305 
0.862442 
0.855326 
0.847611 
0.839222 
0.830082 
0.820129 
0.809309 
0.797485 
0.784694 
0.770972 
0.756151 
0.740230 
0.723140 
0.704991 
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81 


0.0 


0.0 


0.002523 


0.166465 


0.434202 


0.685681 


82 


0. 0 


0.0 


0.001709 


0.145170 


0.463007 


0.665203 


83 


0.0 


0.0 


0.0011 12 


0.125439 


0.441255 


0.643513 


84 


0.0 


0.0 


0.0006 9 ? 


0.107338 


0.419025 


0.620643 


85 


0. 0 


0.0 


0.000416 


0.090901 


0.396407 


0.596714 


86 


0.0 


0.0 


0.000245 


0.076121 


0.373451 


0. 571711 


87 


0.0 


0.0 


0 . 0001 46 


0 . 05 2 c >43 


0.350085 


0.545390 


88 


0. 0 


0.0 


0.000089 


0.051347 


0.326309 


0.517638 


89 


0. 0 


0.0 


0.0 


0.041315 


0.302192 


0.488493 


90 


0.0 


0.0 


0.0 


0.032734 


0.277755 


0.457873 


91 


0.0 


0.0 


0.0 


0.025489 


0.252999 


0.425604 


92 


0.0 


0.0 


0.0 


0.019426 


0.227658 


0.391003 


93 


0. 0 


0.0 


0.0 


0.014421 


0.201531 


0.353733 


94 


0.0 


0.0 


0.0 


0. 010399 


0.175348 


0. 314507 


95 


0.0 


0.0 


0.0 


0.007269 


0.1 49394 


0.274293 


96 


0. 0 


0.0 


0.0 


0.004913 


0.124323 


0.234013 


97 


0. 0 


0.0 


0.0 


0.003195 


0. 100556 


0.194395 


98 


0.0 


0.0 


0.0 


0.001984 


0.07 83 75 


0.155909 


99 


0.0 


0.0 


0.0 


0.001167 


0.058336 


0. 119676 


1 00 


0. 0 


U. 0 


0.0 


0.000653 


0.041607 


0.088149 
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AT I ON NUMBER 1 

NITH ANGLE = 0.5 

AT I ON LATITUDE = 28.50 

ANDARD ATMOSPHERE TEMPEPATURE ADJUSTMENT 
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